Signatures of quantum integrability and nonintegrability 
in the spectral properties of finite Hamiltonian matrices 
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For a two-spin model which is (classically) integrable on a five-dimensional hypersurface in six- 
dimensional parameter space and for which level degeneracies occur exclusively (with one known 
exception) on four-dimensional manifolds embedded in the integrability hypersurface, we investigate 
the relations between symmetry, integrability, and the assignment of quantum numbers to eigen- 
states. We calculate quantum invariants in the form of expectation values for selected operators and 
monitor their dependence on the Hamiltonian parameters along loops within, without, and across 
the integrability hypersurface in parameter space. We find clear-cut signatures of integrability and 
nonintegrability in the observed traces of quantum invariants evaluated in finite-dimensional invari- 
ant Hilbert subspaces. The results support the notion that quantum integrability depends on the 
existence of action operators as constituent elements of the Hamiltonian. 
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I. INTRODUCTION 

An autonomous classical Hamiltonian system with two 
degrees of freedom, specified by some analytic function 
H(p%, qi]P2i Q2) of canonical coordinates, is either inte- 
grable or nonintegrable - tertium non datur. If a second 
integral of the motion can be found, i.e. an analytic func- 
tion I(p\, q\\Vii Q2) which is functionally independent of 
H and satisfies dl/dt = {H, 1} = 0, the system is proven 
integrable. If chaotic trajectories can be detected in the 
phase flow, the system is demonstrably nonintegrable. 
Although it may happen that neither evidence can be 
ascertained in practice for a given H, one or the other 
status is guaranteed to apply. 

A question of long-standing interest has been whether 
an equally clear-cut classification of systems exists in 
quantum mechanics. Translating the criterion of classical 
integrability into quantum mechanics for systems with 
few degrees of freedom opens up loopholes of ambiguity 
that are not easily closed.EJfl Quantum mechanically, a 
second integral of the motion, i.e. an operator I with 
[H, I] = can always be constructed, fo^jexample, via 
time average of an arbitrary operator A.tM Performing 
the time average in the energy representation eliminates 
all off-diagonal matrix elements of A. Which attributes 
of quantum invariants are most sensitive to the integra- 
bility status of the system? 

Quantum chaos research has identified a catalog of at- 
tributes that distinguish quantized nonintegrable from 
quantized integrable systems .ETlI The most widely stud- 
ied distinctive properties pertain to level statistics. How- 
ever, in the extreme quantum limit of a typical model 
system, where the density of energy levels is low, this 
distinction is blurry at best or altogether unrecognizable. 
Only in the energy range where the level density is high, 
which includes the semiclassical regime, do the contrast- 
ing level spacing distributions come into focus. Other 
indicators of quantum chaos are similarly ambiguous. 



One unequivocal discriminant between quantized inte- 
grable and nonintegrable systems was recently identified 
in a study of level crossing.. manifolds in the parameter 
space of a two-spin model.u The system is specified by 
the quadratic Hamiltonian 
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for two quantum spins Si, S2 of equal length y/ 0(0 + 1) 
(a = i, 1, |, . . .). In the classical limit h — > 0, a — > 00, 

%\J a (a + 1) = s, the operators Si turn into 3-component 
vectors, Si = s (sm. -di cos ifi, sin di sin ipi , costfj), and 
Eq. (|l|) then describes the energy function of an au- 
tonomous Hamiltonian system with two degrees of free- 
dom and canonical coordinates pi = s cos i?i , qi = 
<Pi, i = 1,2. The classical integrability condition was 
shown to have the forrnEl 
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Quantum mechanically, the Hamiltonian (|l]) is express- 
ible as a real symmetric block-diagonal matrix, where 
each of the infinitely many finite-dimensional blocks is as- 
sociated with one spin-cr realization of an irreducible rep- 
resentation of the underlying (discrete) symmetry group 
(see Appendix ^) . 

The main conclusions of the level crossing study for 
this system may be summarized as followsn (i) In the six- 
dimensional (6D) parameter space o£-dlJ), level degenera- 
cies occur on smooth 4D structures.ll3(ii) For an invari- 
ant block of H with K levels, this 4D structure consists 
of K — 1 sheets, each representing one twofold [k, k+1] 
level degeneracy in the sequence E% < E% < . . , < Ek- 
(iii) All 4D level crossing sheets are completely embed- 
ded in the 5D integrability hypersurface. (iv) Under mild 
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assumptions, the integrability condition (g) can be deter- 
mined analytically from the conditions of level degener- 
acy in low-dimensional invariant Hilbert subspaces of H. 

These results strongly suggest that the notion of in- 
tegrability remains meaningful for quantum systems de- 
scribed by finite Hamiltonian matrices, notwithstanding 
the fact that there exist universal algorithms for the di- 
agonalization of finite symmetric matrices. 

For a deeper understanding of this subtle notion of 
quantum integrability, we note that classical integrabil- 
ity guarantees the existence of a canonical transforma- 
tion (pi, qi\p2, 92) — > (Ji, J2, #2) to action-angle coor- 
dinates. It transforms the Hamiltonian H(pi,qi;p2,q2) 
and the second integral of the motion I(pi, qi;P2, 92) into 
analytic functions Hc(Ji,J2), Ic(Ji,J2)- Each point 
(Ji, J2) on the action plane specifies a torus in phase 
space. In the nonintegrable case, the actions Ji, J2 are 
only defined for the surviving tori. Since the tori are no 
longer dense anywhere in phase space, no smooth func- 
tions Hc,Ic on Jii J2 exisi-anymore. 

In a companion paperU we have postulated that 
the underlying cause for the embedment of (dj — 1)- 
dimensional level crossing manifolds in a e?/-dimensional 
integrability manifold of the parameter space (with di- 
mensionality d > di) is linked to the existence of ac- 
tion operators as constituent elements of the Hamilto- 
nian. In that study we have demonstrated for two dis- 
tinct model systems the explicit functional dependence 
Hq(Ji, J2), Iq(Ji, J2) of the Hamiltonian and the sec- 
ond integral of the motion on two action operators, and 
compared it to the similar yet different functional depen- 
dence Hc(Ji, J2), Ic(Ji,J2) of the corresponding classi- 
cal invariants on the classical action coordinates. A direct 
comparison was facilitated by the fact that for the model 
systems considered there we knew not only the second 
integral of the motion but also found a set of separable 
canonical coordinates for the description of the classical 
time evolution. 



II. METHOD 

A more indirect but no less compelling method for 
demonstrating the existence of action operators as con- 
stituent elements of the quantum invariants H , / in some 
regions of parameter space, namely on the integrability 
hypersurface, and their nonexistence elsewhere is pur- 
sued here for the two-spin model (|l]). We investigate 
the functional dependence of the eigenvalues of quantum 
invariants on the Hamiltonian parameters, in particular 
across lines demarcating changes in symmetry and/or in- 
tegrability status. 

On the integrability hypersurface (|^), the natural 
quantum numbers of the eigenstates within any invari- 
ant Hilbert subspace of H are the integer pairs (mi, m^) 
specifying the eigenvalues (in units of K) of the action 
operators J\ , J2 ■ Henceforth we call them action quan- 



tum numbers. Elsewhere in parameter space, where 
level crossings between eigenstates of the same param- 
eter space are prohibited, the natural quantum number 
is a single integer, the energy sorting quantum number n. 
What consequences do these conflicting assignments of 
quantum numbers in the two regions of parameter space 
have for the functional dependence of quantum invariants 
on the Hamiltonian parameters? 

Consider the case of a if-dimensional invariant sub- 
space of (0) spanned by the basis given in Appendix [A|. 
The K eigenstates \k),k — 1, . . . , K then form a star of 
orthonormal vectors pointing in oblique directions with 
respect to the coordinate axes. A tiny change of the 
parameters J a , A a causes the star of eigenvectors to ro- 
tate slightly. By monitoring the inner product between 
eigenvectors before and after every infinitesimal param- 
eter change, we can keep track of all eigenvectors along 
the entire loop in parameter space. 

At the same time, we monitor the effect of the gradu- 
ally transforming eigenvectors on the eigenvalues of two 
quantum invariants. For this purpose we choose the en- 
ergy expectation value Ek = (k\H\k) and the expecta- 
tion valuje Ik = (k\A\k), where A is some function of 
the Sf £3 When the Hamiltonian parameters J a , A a are 
varied along a path in 6D parameter space, the vector 
|fc) traces a path on the surface of a X-dimensional unit 
sphere, and the point (Ek, Ik) leaves a trace in the plane 
of invariants. 

What if two eigenvectors are accidentally degenerate 
(Ek = Ek'), which happens when their energy eigenval- 
ues cross each other at some point on the path in param- 
eter space? Generically, the eigenvalues of the second 
invariant are different at the point of level degeneracy 
(Ik 7^ Ik')- We can always choose the second invariant 
such that this is the case. At the crossing point the orien- 
tation of the two eigenvectors is not fixed. However that 
ambiguity is removed if we impose the condition that the 
path of every point (Ek,Ik) in the plane of invariants 
must be continuous. 

We shall see that varying J Q , A a along a closed path 
in parameter space does not guarantee that the trace of 
every eigenstate in the (Ek, ifc)-plane is also closed. It 
may happen, for example, that two eigenvectors trans- 
form into each other in the course of one parameter-space 
loop, thus leaving an open trace in the plane of invari- 
ants, which will be closed only after a second traversal of 
the loop. The two kinds of quantum numbers assigned 
to eigenstates in different regions of parameter space as 
discussed previously, suggest the following scenario. 

(i) If the closed path in parameter space lies entirely 
on the integrability hypersurface, then the traces of all 
eigenstates in the plane of invariants will be closed. Along 
the loop, level crossings occur frequently, but the labeling 
of all eigenstates by the action quantum numbers mi , TO2 
remains valid on every stretch of it. 

(ii) If the path in parameter space lies entirely off the 
integrability hypersurface, the traces of all eigenstates 
will again be closed but for a different reason. Level 
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FIG. 1. Reduced parameter space (J y ,J z ,A) projected 
onto the (J y , J z )-plane. The two dashed lines mark the inter- 
section, 2jf — Jy = 1, of the integrability hyperboloid with 
the plane A = 0. In the integrability plane A — 0, level 
degeneracies of H A1A occur along the dot-dashed lines and 
multiple degeneracies at the symmetry points | J„| = \J Z \ = 1 
marked by the four pentagons. The solid circles represent 
projections of paths with radii \J Jy + Jf = | , \/2, | along 
which we track the quantum invariants Et,Ik- 

crossings are prohibited in this region. All states are 
labeled by the energy sorting quantum number n. That 
label is valid along the entire loop. 

(iii) If the closed path in parameter space consists of a 
leg A on and a leg B off the integrability hypersurface, 
then the conflicting assignment of quantum numbers has 
the consequence that some of the traces in the plane of 
invariants remain open. An eigenstate \k) may undergo 
one or several level crossings on leg A of the path and thus 
end up at a different position in the energy-level sequence 
at the beginning of leg B when the energy-sorting quan- 
tum number kicks in. As the parameters are varied along 
leg B back to their starting values, the point (Ek, Ik) is 
prevented from finding its way back to the original posi- 
tion in the plane of invariants because level crossings are 
now prohibited. 

Not surprisingly, physical reality turns out to be more 
complicated. However, the observations made by this 
method of analysis prove to be highly illuminating in 
regard to the relations between symmetry, integrability, 
and the assignment of quantum numbers. 

III. RESULTS 

To facilitate comparison with results obtained previ- 
ously, we use the same reduced 3D parameter space as 
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FIG. 2. Reduced parameter space (J y ,J z ,A) projected 
onto the (J z , yl)-plane. The solid circles represent paths with 
Jf + A 2 = 0.3712, 0.58 at J y = 0.4 along which we track the 
quantum invariants Ek, Ik- The larger circle is located on the 
integrability hyperboloid. The positions of level crossings of 
H\\a states along that path are indicated by squares. The 
dashed line marks the integrability plane ^4 = 0. 

in Rcf. H It is spanned by—^, J z , A x — A y = 2 A at 
J x = 1, A x + A y = 0, A z = 0£3 The integrability condi- 
tion (H), which becomes 

A(l + J 2 y - 2J 2 Z - 2A 2 ) = 0, (3) 

is satisfied on a 2D surface consisting of the plane A = 
and a hyperboloid with axis at A — 0, J z = 0. Em- 
bedded in this integrability surface are ID level crossing 
manifolds in patterns whose complexity increases with 
the number of levelsrin the invariant (Hilbert) subspaces 
under consideration.tJ 

Individual eigenstates \k) will now be tracked along 
closed paths in this reduced parameter space. Each path 
selected displays distinct characteristic features in the 
traces on the plane of invariants (Ek,Ik)- Here we use 
h = {k\(Sf + ) 2 |fc). We consider invariant (Hilbert) 
subspaces of symmetry class AlA with K = 6, 10 levels 
corresponding to spin quantum numbers a = 4, 5, respec- 
tively (see Appendix |A|) . 

Figure]!] depicts the reduced parameter space projected 
onto the integrability plane A = 0. The dot-dashed lines 
represent the level crossing manifold of H\ lA with K = 
10 levels in the plane A = 0. None of the intersection 
points of two dot-dashed lines involves triple or quadruple 
degeneracies. Each level crossing line can thus be labeled 
[fc, k + 1] by the positions in the level sequence E\ < Ea < 
■ • • < Ek of the two levels involved in the crossing. Ej 

The integrability hyperboloid intersects the integrabil- 
ity plane along the two dashed lines. There exist 30 H\ lA 
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FIG. 3. Energy eigenvalues E^,k — 1, . . . , 10 in the in- 



variant subspace of H A1A as denned in Eq. (A4) and plotted 



ve rsus the angular distance a on the circular path with radius 
yjf + Jl = 0.5 in the plane A = of the reduced parameter 
space (J y , J Z ,A). 




FIG. 4. Closed traces in the (Ek, ifc)-plane of two H A1A 
levels along the circular path with radius -J Jy + Jf = i in 
the plane A = of the reduced parameter space (J y , J Z ,A). 
The traces start at the open squares (a = 0°) in the directions 
indicated. The asterisks mark level crossing points 



level crossing lines on the hyperboloid. These lines inter- 
sect the plane A = at seven points on each dashed line, 
namely on the intersection points with dot-dashed lines 
and on the symmetry points at \J V \ — \J Z \ = 1. The 
solid circles represent projections of paths along which 
we track the quantum invariants Ek, Ik- 

A different projection of the reduced parameter space 
is shown in Fig. |. The larger circle represents a path 
along the intersection of the integrability hyperboloid 
with the plane J y = 0.4. The squares on that circle 
mark the locations where the ten level crossing lines on 
the hyperboloid for H\ 1A intersect the plane J y = 0.4. 
The smaller (concentric) circle represents a path that is 
located in the nonintegrable region of parameter space 
except for the two points where it intersects the integra- 
bility plane A — (dashed line). 



A. Hallmark of integrability 



The first path considered is the circle Jy + J z = \ in 
the plane A = as shown in Fig. [|. This path does not 
come close to any of the symmetry points (pentagons). 
In Fig. H we have plotted the ten levels of H A1A versus 
angular distance a on the circular path. We observe 20 
pairwise crossings between six levels at the angles where 
the path intersects the dot-dashed lines in Fig. 

No instances of level repulsion can be discerned in this 
plot, which is not to say that the a-dependence of ad- 
jacent levels is uncorrelated. Take the six levels near 
the center of the spectrum. They can be divided into 
two groups of three levels undergoing similar oscillations 
along the path. The synchronicity of these oscillations 
is, in fact, a consequence of the (postulated) smooth de- 
pendence of the functions Hq(J\, J2) and Iq(J\, Ja)-.on 
a for this path embedded in the integrability plane.tll 

In Fig. [| we show the traces in the (Ek, ifc)-plane of 



the two eigenstates whose levels undergo four [7,8] cross- 
ings along the path (thick lines in Fig. |J). The traces 
are continuous, closed, and smooth. The square and the 
arrow indicate the starting point and the direction of the 
trace. Every level crossing is represented by two verti- 
cally displaced asterisks, one on each trace. 

It is important to note that the traces remain perfectly 
smooth at the points of level crossing. The level crossings 
have no impact on the eigenvectors, or on the expectation 
values Ik- Every eigenvector loops around and returns 
to its original orientation in Hilbert space. Its path is 
largely unaffected by the presence of other eigenvectors 
which become instantaneously degenerate with it. It is as 
if vectors undergoing level crossings belonged to different 
invariant subspaces. 

The behavior of energy levels as observed in Fig. |^ and 
the properties of traces as seen in Fig. |J reflect what we 
expect for a typical situation in an integrable system with 
two degrees of freedom. The two invariants Ek , Ik are 
functions of two quantized actions J\ , J2 with a smooth 
dependence on the Hamiltonian parameters. The discrete 
values of the actions define the natural quantum numbers 
of all levels, and each eigenstate maintains its identity 
along any path in parameter space notwithstanding the 
presence of level crossings. All traces produced along 
closed paths are therefore closed as well. 

There are two sources of complication forcing on us a 
refinement of this description without undermining the 
postulated link between quantum integrability and ac- 
tion operators. These two complications will be discussed 
next before we investigate the effects of nonintegrability. 



B. Level repulsion near symmetry points 



The second path considered is the circle Jy 
in the integrability plane A = (see Fig. n 
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FIG. 5. Energy eigenvalues E^,k = 1, . . . , 10 in the invari- 
ant subspace of H A1A plotted versus a on the path with radius 
\J Jy + Jz = 1.5 in the plane A = of (J y , J z , A) space. 

makes it different from the previous path is that it passes 
close to the four points \J y \ = \J Z \ = 1, where additional 
degeneracies occur, caused by a higher symmetry. 

The ten levels of H\ 1A versus a are plotted in Fig. ||. 
As in Fig. || for the previous path, we observe 20 level 
crossings, each one associated with a point where the 
circular path intersects one of the dot-dashed lines in 
Fig. [I]. In addition to these crossings we observe instances 
of level collisions at a = mr/2,n — 1,3,5,7, i.e. in the 
vicinity of the symmetry points. 

It is instructive to compare the effects of level crossings 
and level collisions on the traces in the plane of invariants. 
In Fig. H we show again the trace of the point (Ek, Ik) for 
two states that are involved in four [7,8] levels crossings 
(thick lines in Fig. ||), now along the second path. These 
traces exhibit features not seen in Fig. [|. 

We again observe that none of the level crossings leaves 
any mark on the traces, implying that the wave functions 
of the two eigenstates are completely unperturbed by the 
instantaneous level degeneracies (see asterisks). On any 
stretch between successive mutual crossings, both lev- 
els collide with one neighboring level, and each collision 
does have a dramatic effect on the traces of the states 
involved in the collision. Level collisions produce precip- 
itous changes in the second invariant Ik near the closest 
encounter of the colliding levels. The rapid variation of 
expectation values signals a strong perturbation of the 
wave functions in a level collision. The presence of this 
characteristic signature of level collisions is as conspicu- 
ous in the traces shown in Fig. |^ as is their absence in 
the traces shown in Fig. fj. 

In what might be called a hard level collision, the two 
states exchange wave functions in a manner like two bil- 
liard balls exchange momenta in a head-on collision. This 
makes it hard to distinguish a hard collision from a cross- 
ing in a plot such as Fig. |^ because of graphical resolution. 
A plot of one invariant versus the other (Fig. ^|) is much 
more sensitive to that distinction. Here a hard level col- 
lision produces a variation in Ik that looks almost like a 
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FIG. 6. Closed traces in the (Ek, ifc)-plane of two H A1A 
levels along the path with radius yJj~+J£ = 1.5 in the plane 
A = of (Jy,J z ,A) space. The traces start at the squares 
(a — 0°) in the directions indicated. The asterisks mark level 
crossing points. 

discontinuity. 

The phenomena observed in Figs. [| and || are not 
in contradiction with the assertion that the invariants 
Ek,Ik are functions of two quantum actions. It tells 
us, however, that the dependence of these functions on 
the Hamiltonian parameters is singular at the symmetry 
points of H . The phenomenon of level repulsion in the 
immediate vicinity of symmetry points is then caused by 
invariants pertaining to the higher symmetry and by the 
associated additional level degeneracies. 

The traces of all levels depicted in Fig. la are closed as 
were all traces of the levels shown in Fig. |3j The implica- 
tion is that the number of crossings between any pair of 
levels must be an even number. The fact is that neither 
the level crossings nor the level collisions can cause any 
confusion in the labeling of the levels by action quantum 
numbers along a path in the integrability plane A = as 
long as it avoids the points \J y \ = \J Z \ — 1 of higher sym- 
metry with symmetry induced level degeneracies. Each 
eigenstate maintains its identity along such paths, or so 
it seems. 



C. Open traces caused by a change in symmetry 

The third path considered is the circle J y + J z = 2 at 
A = (see Fig. |l|). It is embedded in the integrability 
plane and passes through the points \J y \ — \J Z \ = 1. The 
impact of these symmetry points on the energy levels is 
depicted in Fig. 0. What were level collisions in Fig. || 
have now turned into additional level crossings. At the 
symmetry points, the ten levels combine into a singlet, 
a doublet, a triplet, and a quadruplet. No instances of 
level repulsion are observable anymore. 

The absence of level collisions along this path is con- 
firmed by a study of the traces in the (Ek, Ik) -plane. In 
Fig. we show the traces of the two states that again 
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FIG. 7. Energy eigenvalues Ek, k, . . . , 10 in the invariant 
subspace of H A1A plotted versus a on the path with radius 
\J Jy + Jz = V% in the plane A = of (J y , J z , A) space. 



start in the seventh and eighth positions of the level se- 
quence. Gone are the rapid near- vertical displacements 
which we have identified in Fig. |^ and which were caused 
by level collisions. The traces in Fig. || are as unaffected 
by the new symmetry-induced level crossings as they are 
oblivious of crossings elsewhere in the integrability plane. 

However, a striking new feature makes its appearance 
in Fig. ||. The traces do not close in themselves after one 
loop around the circular path in parameter space. The 
eighth level becomes the seventh level after one loop, and 
then turns into the second level after two loops. Only 
after the third loop does it end up in the original eighth 
position of the level sequence. 

In Fig. [?] the three levels involved in that loop are 
drawn as thick lines. Inspection shows that there are 
two further groups of three states which transform into 
each other as the parameter values loop around the circle. 
That leaves one state (near the center of the spectrum) 
whose trace closes in itself after one loop. 

Is this phenomenon of levels transforming into each 
other compatible with the notion that the invariants are 
functions of the quantized actions with a smooth depen- 
dence on the Hamiltonian parameters? Yes if we al- 
low the dependence on the parameters to be singular at 
points of higher symmetry within the integrability mani- 
fold. The presence of such singularities was already sug- 
gested by the level collisions observed in Figs. [| and g. 
The results of Figs. and || confirm the singular param- 
eter dependence from a different vantage point. 

When we start with the second path in parameter 
space (Sec. IIIB) and increase the radius of the circle 



gradually toward that of the third path, we observe a 
gradual hardening of the level collisions near the sym- 
metry points. The hardening is characterized by increas- 
ingly sharp curvatures in the graphs of Ek versus a (Fig. 
^) and by increasingly rapid vertical variations in the 
graphs I k versus E k (Fig. g). 

In the limiting case of this path, the sharply curved but 
smooth bends in the graph E k versus a turn into cusps, 




FIG. 8. Open traces in the (Ek, 7fc)-plane of three H A1A 
levels along the path with radius Jy + Jf = \/2 in the plane 
A — of (J y , Jz,A) space. The traces start at the squares 
(a — 0°) in the directions indicated The asterisks mark level 
crossing points. 



and the fast but smooth vertical variations in the graphs 
Jfe versus Ek turn into discontinuities. An infinitely hard 
level collision is indistinguishable from a level crossing. 
In Figs. |^ and || smooth segments of graphs between 
singularities that belong to different colliding levels are 
rejoined to form entirely smooth graphs of crossing levels. 

Hence, if we insist that all levels maintain their identity 
along any closed path in the integrability plane A = 0, 
we must interpret all level crossings that take place at the 
points of higher symmetry, \J y \ = \ J Z \ = 1 , as infinitely 
hard level collisions. All the evidence accumulated thus 
far still supports the existence of the functions Hq(Ji, J2) 
and Iq(Ji, J2) with a smooth parameter dependence on 
the integrability manifold, provided we allow for singu- 
larities at points of higher symmetry. 

Before we discuss the strongly contrasting properties 
of quantum invariants along paths that are not fully em- 
bedded in the integrability manifold of (|l|), we should 
report on yet another feature that complicates the inter- 
pretation of the integrable cases. 



D. Open traces caused by topology 

The circle A 2 + J 2 = 0.58 with center at J y = 0.4 
is the fourth path along which we study the behavior 
of quantum invariants. This path represents a circular 
section of the integrability hyperboloid (||) (see Fig. ||). 
Like the first path considered, it does not pass near any 
point in parameter space where symmetry induced level 
degeneracies occur. 

The angular dependence of the six H\ lA levels, de- 
picted in Fig. |9[ does indeed not show any level collisions 
just as was the case in Fig. ||for the first path. All levels 
undergo several crossings along this path, and none of 
the crossings has any noticeable effect on the quantum 
invariants Ek,Ik plotted in Fig. 10. 
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FIG. 9. Energy eigenvalues E k ,k = 1, . . . , 6 in the invariant 
subspace of H A1A plotted versus a on the path with radius 
\f A 2 + Jf = V0.58 at J y — 0.4 embedded in the integrability 
hyperboloid of {J y , J z , A) space. 
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FIG. 10. Open traces in the (E k ,I k 



-plane of all six H A1A 



levels along the path with radius A 2 + J 2 



V0.58 at 
A) space. 

The traces start at the squares (a = 0°) in the directions 
indicated. The asterisks mark level crossing points. 



J y = 0.4 on the integrability hyperboloid in (J y , J, 



FIG. 12. Closed traces in the (E^, 7fc)-plane of all six H A1A 
levels along the path with radius V A 2 + J 2 = \/0.3712 at 
J y — 0.4 off the integrability hyperboloid in (J y , J z , A) space. 
The traces start at the the squares (a = 0°) in the directions 
indicated. 



Nevertheless, there is a major difference between the 
evolution of eigenstates along these two paths. Each one 
of the six levels shown in Fig. ^ transforms into a different 
level in the course of one loop of the path around the 
integrability hyperboloid. It takes three loops for every 
eigenstate to return to its original position in the level 
sequence. On the plane of invariants this phenomenon is 
reflected in open traces that connect to form two rings of 
three segments each as shown in Fig. [l(]. The two sets 
of levels are distinguished by line thickn ess. 

Unlike in the previous situation (Sec. Ill C| ) , here the 
open-trace phenomenon cannot be attributed to a change 
of symmetry along the path. What distinguishes the first 
path, where open traces do not occur from the fourth 
path, where they do occur, is that only the former can be 
shrunk to a point without leaving the integrability mani- 
fold. Hence the multiple connectedness of the integrabil- 
ity hyperboloid forces us to allow for functions Hq ( Ji , J 2 ) 



and Iq ( Ji , J2) whose dependence on the Hamiltonian pa- 
rameters is still smooth but multiple-valued. 

With these concessions, the signature properties of 
quantum integrability postulated above remain fully in- 
tact. The quantum invariants E k ,I k exhibit strongly 
contrasting features when observed along paths that are 
not embedded in the integrability manifold. Visualizing 
these differences does not depend on a statistical anal- 
ysis. They are unmistakenly identifiable in systems for 
systems with very few levels. 



E. Level repulsion due to nonintegrability 

For a direct comparison with the previous situation, we 
now choose a circle with the same center as the fourth 
path and a somewhat smaller radius, J 2 + A 2 = 0.3712. 
This fifth path lies off the integrability manifold except 
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for two points where it intersects the integrability plane 
A = (see Fig. ||). However, no level degeneracies occur 
at these intersection points. 

The six H\ 4 levels versus a along the fifth path are 
plotted in Fig. flU Even though the resulting pattern is 
vaguely similar to that observed in Fig. the differences 
are clear-cut. All level crossings have turned into level 
collisions. 

Most of the collisions are fairly soft. The two hard- 
est collisions are barely resolved as such on the scale of 
Fig. |l2]. None of the levels transform into each other 
any more. The levels are now naturally labeled by the 
energy sorting quantum number. Each open segment of 
the traces shown in Fig. |l0| has turned into a closed 
trace. All level collisions, especially the hard ones, leave 
the characteristic marks on the traces in the form of a 
rapidly varying second invariant Ik- 

If we were to move the fifth path closer to the integra- 
bility hyperboloid by increasing its radius (see Fig. ||) , we 
could observe a gradual hardening of all level collisions. 
The level configurations as shown in Fig. [fl] would in- 
creasingly resemble those in Fig. || The traces as shown 
in Fig. [l^, however, would remain very different from 
those pertaining to the integrable case (Fig. |lCi| ). 

Only in the limiting case where the fifth path merges 
with the fourth path would the closed traces of the nonin- 
tegrable model break into segments connected by vertical 
lines. The ends of each segment would then rejoin ends of 
other segments to form the smooth rings of open traces 
shown in Fig. [l^. 

Similar observations are made upon lifting the first 
path off the integrability plane A = to a plane at A 7^ 0. 
All the level crossings that exist in Fig. [|, for example, 
turn into level collisions. The closed traces such as shown 
in Fig. ^ break into pieces whose ends rejoin via near ver- 
tical lines into a new set of closed traces. 

Along the second path we had observed (in Fig. [|) level 
crossings (due to integrability) and level collisions (due to 
nearby points of higher symmetry) . Lifting this path off 
the integrability plane again removes all level crossings 
and results in a set of closed traces. The characteristic 
marks of level collisions on the traces in the (Ek,Ik)~ 
plane are the same no matter whether they are caused 
by a reduced symmetry or by nonintegrability. 

Lifting the third path off the integrability plane has 
the same effects on the level crossings attributed to inte- 
grability and the level crossings attributed to the higher 
symmetry at selected points in parameter space (Fig. [7]). 
All are removed indiscriminately. 



F. Open traces caused by nonintegrability 

The conflicting assignments of quantum numbers to 
eigenstates for parameter values on and off the integra- 
bility manifold is most compellingly documented when 
we pick a path in parameter space that is only partially 




360 



FIG. 13. Energy eigenvalues Ek,k — 1, . . . , 10 in the in- 
variant subspace of H A1A plotted versus a on the path with 
radius y/Jf+J% = 0.5 at A = 0.3 cos 2 (a/2) in (J y ,J z ,A) 
space. 




FIG. 14. Open traces in the (Eh, Ifc)-plane of the two H A1A 
levels which undergoes one [5,6] crossings along the path with 
radius a/ J 2 + Jf = 0.5 at A = 0.3 cos 2 (a/2) in (J y ,J z ,A) 
space. The traces start at the squares in the direction indi- 
cated. The asterisk on each trace marks the level crossing 
point at a = 180° . 



embedded in the integrability manifold. 

The sixth path considered in this study of quantum 
invariants is a modification of the first path (Sec. [II A) 
with the same projection in Fig. [j]. Whereas the first 
path was embedded in the integrability plane A = 0, the 
sixth path has a variable height relative to that plane: 
A(a) — 0.3 cos 2 (a/2). It touches down to the integrabil- 
ity plane at a single point (a = 180°), where a [5,6] level 
crossing takes place. 

Along this path there exist no other level crossings. 
All the other crossings that existed in Fig. || for the first 
path are now replaced by level collisions (see Fig. |l3| ). 

The inevitable consequence of having a single level 
crossing along a closed path in parameter space is the 
existence of a pair of open traces in the plane of invari- 
ants, namely the traces of the states that undergo the 
[5,6] crossing at a = 180°. These traces are shown in 
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Fig. The ends of the solid and dashed lines form a 
single loop, which is traced in the direction indicated. 

What causes here an open trace in the plane of invari- 
ants is obviously akin to what had cause d an open trace 
in the situation described in Sec. Ill C . In both cases 



two levels cross once due to particular circumstances at 
one point of the path, and are thus prevented from cross- 
ing back to their original position in the level sequence 
on the remaining stretch of the path. In Sec. Ill C the 



particular circumstance was a higher symmetry, here it 
is integrability. 



IV. INTERPRETATION 

The study of quantum invariants along closed paths 
through parameter space indicates that a change in sym- 
metry and a change in integrability status produce re- 
lated phenomena. In some dynamical systems, the con- 
servation laws that guarantee integrability are direct con- 
sequences (via Noether's theorem) of continuous symme- 
tries. Switching from integrability to nonintegrability is 
then accompanied by a reduction in symmetry. 

In the two-spin model ([l]), the presence of a (contin- 
uous rotational) 0(2) or higher symmetry in spin space 
does indeed imply the existence of a second integral of 
the motion, namely the component of the total spin along 
the symmetry axis, and integrability is guaranteed. How- 
ever, a second integral of the motion was shown to exist 
for certain parameter values even in the absence of a con- 
tinuous rotational symmetryB Does integrability in that 
case indicate the presence of a hidden symmetry? 

Classical integrability guarantees that the Hamilto- 
nian (|l|) can be expressed as a function of the two ac- 
tion variables: H — Hc(Ji, J^)- The cyclical nature of 
the angle coordinates thus implies that He is invariant 
with respect to continuous rotation-like transformations 
in phase-space. Since this is not related to a continuous 
symmetry in configuration space, it is appropriate to call 
it a hidden symmetry. 

For a description of the impact of symmetries on the 
level spectrum of the quantum two-spin model, it is useful 
to distinguish three kinds of symmetry: discrete symme- 
tries, continuous symmetries, and hidden symmetries. 

Discrete symmetries have no bearing on the classical 
integrability property, but they do affect the shapes of 
phase-space trajectories. Quantum mechanically, they 
divide the Hilbert space of H into invariant subspaces. In 
general, this does not result in symmetry-induced level- 
degeneracies, but it does lead to accidental degeneracies 
between levels belonging to different invariant subspaces. 
Such level crossings exist independently of whether or not 
H is integrable. 

Hidden symmetries, which guarantee classical inte- 
grability, cause additional accidental level degeneracies, 
namely between states within one of the invariant sub- 
spaces pertaining to any existing discrete symmetry. 



Continuous symmetries, in essence, combine the effects 
of the discrete and hidden symmetries, and allow acciden- 
tal inter-subspace degeneracies. In addition to these ef- 
fects, continuous symmetries (sometimes in tandem with 
discrete symmetries) produce level degeneracies of a per- 
manent nature, the so-called symmetry-induced level de- 
generacies. 

There exists a hierarchy of symmetries in the two-spin 
model ([!]) : (SO) In the absence of any symmetry, there are 
no level degeneracies. All levels will collide when Hamil- 
tonian parameters are varied. This situation can be re- 
alized by an external magnetic field. (SI) The existence 
of discrete symmetries alone produces finite-D invariant 
Hilbert subspaces. Level crossings exist between states 
belonging to different subspaces. Levels within any sub- 
space collide. (S2) The existence of hidden symmetries in 
addition to discrete symmetries produces level crossings 
between states in the same invariant subspace. (S3) The 
continuous symmetries produce permanent degeneracies 
in certain regions of parameter space. 

There exists a hierarchy of level collisions which corre- 
sponds to the hierarchy of symmetries. (SI— >S0) Inter- 
subspace level crossings in the presence of discrete sym- 
metries turn into level collisions when discrete symme- 
tries are removed. (S2— >1) Intra-subspace level crossings 
turn into level collisions when the hidden symmetries 
are removed, i.e. when the integrability is destroyed. 
(S3— >S2) Symmetry-induced level degeneracies associ- 
ated with a continuous symmetry are removed outside 
the range of that symmetry irrespective of the presence 
or absence of the hidden symmetry. 

Some level crossings along paths through symmetry 
points in parameter space turn into level collisions along 
nearby paths that miss the symmetry point. Other level 
crossings are insensitive to whether the path hits or 
misses the symmetry point. They are the product of 
the hidden symmetry. 

All phenomena observed in the quantum invariants 
Ek,Ik along closed paths on, off, and across the inte- 
grability manifold, indicate that the effects of a change 
in integrability status are akin to the effects of a change 
in symmetry. All observations point to the existence of 
a hidden symmetry that accompanies quantum integra- 
bility. 

In the classical limit, this hidden symmetry manifests 
itself in phase space when viewed from a particular coor- 
dinate system - the action-angle coordinates. The same 
hidden symmetry must also exist in the quantum system, 
but only on the integrability manifold. Even though non- 
integrability is not to be taken literally in the quantum 
case, the presence or absence of that hidden symmetry 
has consequences that are equally clear-cut as in the clas- 
sical limit. 
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APPENDIX A: DISCRETE SYMMETRIES 

The (discrete) symmetry group relevant for the general 
2-spin Hamiltonian (jl]) is Z^CgiS^, where D2 contains the 
three twofold rotations C%, a = x,y, z about the coordi- 
nate axes, and S2 the permutations of the two spins. The 
eight irreducible representations of D2 ® S2 are named 
A1S, A1A, BIS, B1A, B2S, B2A, B3S, B3A, where S 
(A) stand for (anti-)symmetric under permutation and 
Al, Bl, B2, B3 for (Cf,Cf,Cf) = 0.1,1), (1,-1,-1), 
(-1, 1, -1), (-1, -1, 1), respectivelyBB 

TABLE I. Symmetry-adapted basis vectors for integer 
a. The local spin quantum numbers satisfy the relations 
< h < a, —h < h <h- The subspace dimensionality K is 
|(er + l)(<r + 2) for A1S, \o(o - 1) for A1A, and \a(a + 1) 
for the other six classes. 



A1S 


\00) ,(\h,h) + \-k,-h)) /y/2, 
(\h,-h} + \-h,k))/V2, 
(\h,h) + \-h,-h) + |Z 2 , Zi) + |-Z 2 ,-Zi))/2 
Zi + h even. 


A1A 


(|Zi,Z 2 ) + \-h,-k) - \h,h) - \-h,-h))/2 
Zi + h even. 


BIS 


(|Zi,Zi) - |-/i,-/i)) /V2 
(|Ii,Ia> - Hi, -fa) + |fa,Zi) - |-fa,-Zi))/2 
h + h even. 


B1A 


(|li,-/i> - \-h,h))/V2 
(|Zi,Z 2 ) - |-/i,-/2> - |fa,Zi) + |-/2,-/i)) /2 
Zi + h even. 


B2S 


(|Ii,Ia> - 1-/1,-^2) + |I a ,Ii> - 1-/2,-/1)) /2 
/1 + I2 odd. 


B2A 


(|/i,/ 2 ) - 1-/1,-/2) - Mi) + |-/2,-/i)) /2 

Zi + h odd. 


B3S 


(|/i,/ 2 ) + 1-/1,-/2) + |/ 2 , /1) + |-/2,-/i)) /2 
Zi + Z2 odd. 


B3A 


(|/i,/ 2 ) + |-/i,-/2) - |/ 2 , /1) - |-/2,-/i)) /2 
Zi + h odd. 



The basis vectors with transformation properties cor- 
responding to the eight different irreducible representa- 
tions R are listed in Table | for integer a and in Table @ 
for half-integer a. The Hamiltonian matrix can then be 
expressed in the form 



H 



(Al) 



with blocks of dimensionalities K = 1, 3, 6, 10, . . . in 16 
different realizations, two for each symmet:™ class (one 
with integer a and one with half- integer a)x3 



TABLE II. Symmetry-adapted basis vectors 
for half- integer a. The local spin quantum numbers satisfy 
the relations < Zi < a, — Zi < Z 2 < Zi. The subspace dimen- 
sionality if is I (4<r 2 — 1) for symmetric and | (2<r + 1) (2a + 3) 
for antisymmetric representations. 



A1S 


(\U U) - 1— Zi -U) + l/o U) — \—U —U\\ 12 
^1+^2 even. 


A1A 


(|/i,-/i)-|-/i,/i))/V2 
Cl/i ;,\ — I—/, — ;,\ — 1/., /-,\ 4- 1—/, — ],)) 12 

\\ b 1 1 b 2 / i>1 ; i>2/ \ b 2 ) b 1 / \ \ b 2i t>l / J / ^ 

h + li even. 


BIS 


n/i — /1 \ a- i— /1 /1 \^ up) 

(|/i, h) + 1-/1,-/2) + |/ 2 , Zi) + I-/2, -Zi)) 12 
h + h even. 


B1A 


CIZ1 lo) + l-/-i -M — I/9 h\ — I— U — Ml/2 
h + I2 even. 


B2S 


(|/i,/i) + \-h,-h})/s/2 
(|/i,/ 2 ) + |-/i,-/2) + |/ 2 , Ji) + I-/2, 12 
li + h odd. 


B2A 


(|/i,/ 2 ) + 1-/1,-/2) - |/ 2 , /1) - |-/ 2 ,-/i))/2 
Zi + h odd. 


B3S 


(|/i,/i) - |-/i,-/i))/V2 
(|/i,/ 2 ) - |-/i,-/2) + |*2, Zi) - |-/2,-/i» /2 
h + h odd. 


B3A 


(|/i,/ 2 ) - 1-/1,-/2) - |/ 2 , h) + |-Z 2 ,-/i))/2 
Zi + h odd. 
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